A general solution in the Newtonian limit of f(R) - gravity 
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We show that any analytic /(i?)-gravity model, in the metric approach, presents a weak field 
limit where the standard Newtonian potential is corrected by a Yukawa-like term. This general 
result has never been pointed out but often derived for some particular theories. This means that 
only f(R) — R allows to recover the standard Newton potential while this is not the case for other 
relativistic theories of gravity. Some considerations on the physical consequences of such a general 
' solution are addressed. 
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c : 

/(i?)-gravity can be considered a reliable mechanism to explain the cosmic acceleration by extending the geometric 
sector of field equations [l|, 0, H, 0, M, B 0> H, B EH EH EH EH EH EH ■ There are several physical and mathematical 
IO I motivations to enlarge General Relativity (GR) by these theories. For comprehensive reviews, see [H, E3 EH ■ 

Besides, dealing with such extended gravity models at astrophysical scales, one faces the emergence of corrected 
' gravitational potentials with respect to the Newton one coming out from GR. This result is well known since a long 
Q" 1 time [IH, [2(| but recently it has been pursued to carry out the possibility to explain the flatness of spiral galaxies 

rotation curves and the potential of galaxy clusters without adding hu ge a mounts of dark matter [22|, HH . 
(3JT)! Other issues as, for example, the observed Pioneer anomaly problem [24 [25j can be framed into the same approach 
[2^ | and then, apart the cosmological dynamics, a systematic analysis of such theories urges at short scales. 

In this letter, we discuss the Newtonian limit of analytic /(i?)-gravity models deriving a general solution for the 
gravitational potential. 

The discussion about the short scale behavior of higher order gravity has been quite vivacious in the last years since 
GR shows its best predictions just at the Solar System level. As matter of fact, measurements coming from weak 
field limit tests like the bending of light, the perihelion shift of planets and the frame dragging experiments represent 
inescapable tests for any theory of gravity. Actually, in our opinion, there are sufficient theoretical predictions to state 
t-H ■ that higher order theories of gravity can be compatible with Newtonian and post-Newtonian prescriptions [28| since 
the standard Solar System tests can be evaded by several classes of them [27| ■ 

Nevertheless, up to now, the discussion on the weak field limit of f(R) - theories is far to be definitive and there 
are several papers claiming for opposite results [29|, [3(| [U, HH HH, [3, [35|, HH , or stating that no progress has been 
, reached in the last forty years due to several misconceptions in the various theories of gravity [37l . [38j] . 

In particular, people approached the weak limit issue following different schemes and developing different parame- 
terizations which, in some cases, turn out to be not necessarily correct. 

The present analysis is based on the metric approach, developed in the Jordan frame, assuming that the observations 
are performed in it, without resorting to any conformal transformation (see also |40(). This point of view is adopted 
in order to avoid dangerous variable changes which could compromise the correct physical interpretation of the 
results. As we will see below, a general gravitational potential, with a Yukawa correction, can be achieved in the 
Newtonian limit of any analytic /(i?)-gravity model. From a phenomenological point of view, this correction allows 
to consider as viable this kind of models even at small distances, provided that the Yukawa correction turns out to 
be insignificant in this approximation as in the so called "chameleon mechanism" [39j . From the point of view of 
the post-Newtonian corrections, such a solution implies a redefinition of the coupling constants in order to fulfill the 
experimental prescriptions. 

As matter of fact, if one evaluates what the corrections are to the Newtonian potential coming from a modified 
gravity model in a post-Newtonian regime, it is necessary to take into account corrections both at the second order and 
at the fourth order in the perturbation expansion of the metric. Furthermore, this analysis turns out to be coherent 
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only if a small r regime is adopted. In this limit, in fact, the potential can be retained reasonably Newtonian-like in 
relation with experimental results. 

Let us start from a general fourth order gravity action : 



A = 



f{R) + XL, 



(1) 



where f(R) is an analytic function of Ricci scalar, g is the determinant of the metric c^ t „, X = is the coupling 

constant and C m describes the standard fluid-matter Lagrangian. It is well known that such an action is the straight- 
forward generalization of the Hilbert-Einstein action of GR obtained for f(R) = R. Since we are considering the 
metric approach, field equations are obtained by varying ([!]) with respect to the metric : 



and the trace is 



- \fg»» - f',^ + g MV nf = f Tp, , (2) 



3D/' + /'i?-2/=yT. (3) 



Tau = p= ^_^Zj^!!L i s the energy momentum tensor of matter (T is its trace), /' = ^-r^- and □ = . a ' a is the 
\j— g og^ v dR 

d'Alembert operator. 

Actually, as discussed in (4(j, we deal with the Newtonian and the post-Newtonian limit of f(R) - gravity adopting 
the spherical symmetry. The solution of field equations can be obtained considering the metric : 



ds 2 



g rTT dx cr dx T — goo(x° ,r)dx° 2 + g rr (x° ,r)dr 2 — r 2 d£l 



(4) 



where 



„o 



(•! and dQ is the angular element. 



In order to develop the Newtonian limit, let us consider the perturbed metric with respect to a Minkowskian 
background g^ — rj^v + h„v The metric entries can be developed as: 



fftt (t,r)~l+ 5t ( t 2) (t,r)+ 5t ( t 4) (i,r) 
g rr (t,r) ~ -1 + g { r){t,r) 
9ee(t,r) = -r 2 

. g<t>S,r) 



(5) 



-r 2 sin 2 ( 



where we are assuming c = 1 , x° = ct — > t and applying the formalism developed in [40( • 

Since we want to obtain the most general result, we does not provide any specific form for the /(i?)-Lagrangian. 
We assume, however, analytic Taylor expandable f(R) functions with respect to a certain value R = Rq : 



f(R) = E LA 7^(R - Rn) n * fo + hR + hR 2 + hR 3 + • (6) 

n 

In order to obtain the post-Newtonian approximation, one has to insert expansions ([5]) and ([5]) into field equations 
([2]) - §5$ and expand the system up to the orders O(0), 0(2) e 0(4). This approach provides general results and 
specific (analytic) theories are selected by the coefficients fi in Eq.©. It is worth noting that, at the order 0(0), the 
field equations give the condition /o = and then the solutions at further orders do not depend on this parameter as 
we will show below. If we now consider the 0(2) - order approximation, the equations system ^ in vacuum, results 
to be 
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' hrR (2) 


- 2f l9 % r + 8/ 2 i? ( r 2) - 


hrg%}„ + = 




— %fl9rr,r + 8/2 R^? ~ 


- firgulrr = 


Of « (2) 
ijxg rr 


r[hrR^ hg[% - 


/iS$r + 4/ 2 i?fr } + 4/ 2 rfl$] = 


hrRW 


+ 6/ 2 [2i? (2) + rR {2) r ] = 


= 


s %Qrr ~\~ 


r[2flgi - rRV) + 2g% 


r + rfl^] = 



It is evident that the trace equation (the fourth in the system (J7J), provides a differential equation with respect to 
the Ricci scalar which allows to solve exactly the system (JJJ at 0(2) - order. Finally, one gets the general solution: 



(2) 

9tt 



(2) 
'frr 



5o 



Y 
fir 



Si(t)e- 



3Jr 



<?2(t)e'' v 
6(~«) 3/ 



y 



Jj(2) _ gi(t)e- 



3fr 



<5 2 (t)V=?e 



ggjt)[|rW31e 
6£ 2 r 



(8) 



2£r 



/l 

where £ = 7-7-, /1 and / 2 are the expansion coefficients obtained by Taylor developing the analytic f(R) Lagrangian 
6/ 2 

and Y is an arbitrary integration constant. When we consider the limit / — ► R, in the case of a point-like source of 
mass M, we recover the standard Schwarzschild solution with Y = 2GM. Let us notice that the integration constant 
5q is dimensionlcss, while the two arbitrary functions of time Si(t) and <5 2 (i) have respectively the dimensions of 
lenght~ l and lenght~ 2 ; £ has the dimension lenght~ 2 . The functions of time 6i(t) (i = 1,2) are completely arbitrary 
since the differential equation system (|7|) contains only spatial derivatives and can be settled to constant values. 
Besides, the integration constant 5o can be set to zero since it represents an unessential additive quantity for the 
potential. 

We can now write the general solution of the problem considering the previous expression ([5]). In order to match 
at infinity the Minkowskian prescription for the metric, we can discard the Yukawa growing mode in (JSJ and then we 
have : 



ds 2 = 



1GM _ Si{t)e" 
fir 



3£r 



dt 2 - 



2GM ai(t)(rV=f+l)e 



fir 



3{r 



dr 2 - r 2 dn 



(9) 



I R 



Si(t)e 



At this point, one can provide the solution in term of the gravitational potential. The first of © gives the second 
order solution in term of the metric expansion (see the definition (O). This term coincides with the gravitational 



potential at the Newton order. In particular, since g±t = 1 
of /(i?)-gravity, analytic in the Ricci scalar R, is 



2$ 



grav 



= 1 



(2) 

gu 5 



the general gravitational potential 



grav 



GM 6i(t)t 
+ 



-rV=T 



6£r 



(10) 



This general result means that the standard Newton potential is achieved only in the particular case f(R) = R while it 
is not so for generic analytic f(R) models. Eq. (fTTJ)) deserves some comments. The parameters /i j2 and the function 6% 
represent the deviations with respect the standard Newton potential. To test such theories of gravity inside the Solar 
System, we need to compare such quantities with respect to the current experiments, or, in other words, Solar System 
constraints should be evaded fixing such parameters. On the other hand, such parameters could acquire non-trivial 
values (e.g. /1 7^ 1, Si(t) 7^ 0, £ 7^ 1) at scales different from the Solar System ones. 
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In conclusions, we have derived the Newtonian limit for analytic /(-R)-gravity models. The approach which we 
have developed can be adopted also to calculate correctly the post-Newtonian parameters of such theories without 
any redefinition of the degrees of freedom by some scalar field. This last step could be misleading in the weak field 
limit. In the approach presented here, we do not need any change from the Jordan to the Einstein frame [4ll l42l] . 

Apart the possible shortcomings related to non-correct changes of variables, any f(R) theory can be rewritten 
as a scalar-tensor one or an ideal fluid, as shown in [43|, 0, Ufa]- In those papers, it has been demonstrated that 
such different representations give rise to physically non-equivalent theories and then also the Newtonian and post- 
Newtonian approximations have to be handled very carefully since the results could not be equivalent. In fact, the 
further geometric degrees of freedom of f(R) gravity (with respect to GR), namely the scalar field and/or the further 
ideal fluid, have weak field behaviors strictly depending on the adopted gauge which could not be equivalent or difficult 
to compare. In order to circumvent these possible sources of shortcomings, one should states the frame (Jordan or 
Einstein) at the very beginning and then remain in such a frame along all the calculations up to the final results. 
Adopting this procedure, arbitrary limits and non-compatible results should be avoided. 

As we have seen, the solution relative to the gtt metric component gives the gravitational potential corrected by 
Yukawa-like terms, combined with the Newtonian potential. In relation to the sign of the coefficients entering gtt, 
one can obtain real or complex solutions with different physical meanings. This degeneracy could be removed once 
standard matter is introduced into dynamics [39( . On the other hand, the growing Yukawa mode can be neglected as 
soon as the correct asymptotic Minkowskian limit is required. Such considerations acquire an important meaning as 
soon as f(R) gravity is considered from the viewpoints of stability and Cauchy problem. It is possible to show that 
the correct Newtonian limit of the theory is strictly related to the stability of the ground state solution (4|| |4?l [48|, 149[ 
and this fact determines the range of viable parameters in the potential f) 1 0[) . 
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